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در مطالعه حاضر با استفاده از روش محاسباتی به نام بررس یگروهی داده (۱۱۵۱1:6 1016 ]0 1۱61۵۵ ۲0۵۷) مدل موثر برای پیش‌بینی میزان 
پیشروی واکنش نیتریده‌کردن سیلیسیم به‌روش پیوند واکذشی, تو سعه داده شده است. ۲۱۸۲ داده ا زآزمایش‌های انجام شده در این تحقیّق و اطلاعات 
موجود د رکارهای گذشته استخراج شده است. تأثیرات مهم‌ترین پارامترها مثل زمان, دماء در صد گاز نیتروژن, اندازُ قرص و اندازة ذر سیلیسیم اولیه 
بررسی شده است. واکنش نیتریده‌کردن سیلیسیم در دمای ۱۵۷۳ درجه کلوین انجام شد و نتای جآزمایش با مقادیر حاصل از مدل مقایسه و با استفاده از 
پارامترهای حطای آماری عملکرد مدل با حطای ۰/٩۲‏ تأیید شد. همچنین عملکرد دقیق و موثر مدل 211111 در برابر مدل‌های تحلیلی قبلی مقایسه 
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مقد مه 

سیلیسیم نیترید به‌عنوان یک از مهم‌ترین سرامیک‌های 
مهندسی. کاربرد وسیعی در زمینه مواد دما بالا دارد. مزایای 
زیاد مثل مقاومت بالا در برابر شوک حرارتی و خوردگی؛ 
استحکام بالاء چگالی پایین در مقایسه با ترکیبات فلزی و 
زیست سازگاری باعث شده است که از این ماده در ابزار 
برش. توربین‌های گازی. موتور دیزل. یاتاقان‌ها و 
کامپوزیت‌های زیست سازگارپذیر استفاده شود [1-5]. 
چندین روش مثل پرس گرم سیلیسیم نیترید. پرس گرم 
ایزواستا تیک سیلیسیم نیترید. سیلیسیم نیترید تفت جوش‌شده 
سیلیسیم نیترید تفت‌جوش‌شده در فشار جزئی گاز سیلیسیم 
نیترید پیوند واکنشی و سیلیسیم نیترید پیوند واکنشی 
تفت جوش‌شده برای تولید بدنه‌های سیلیسیم نیتریدی وجود 
رت اه ای ظالخه وشن رسای سطالته اه 
است. 

روش پیوند واکنشی که با عنوان نبتریده‌کردن مستقیم 
نیز شناخته می‌شود. یکی از روش‌های مهم اقتصادی در تولید 
بدنه‌های سیلیسیم نیتریدی است. این روش قطعه را بسیار 
نزدیک به هندسه نهایی تولید می‌کند و نسبت به سایر روش‌ها 
عملیات اضافی بسیار کمی نیاز دارد. همچنین هزین اجرای 
روش پیوند واکنشی رابطه‌ای مستقیم با زمان تولید قطعه دارد. 
در این روش سیلیسیم نیترید از طریق واکنش نیتروژن و 
سیلیسیم در دمای بین ۱۲۵۰ درجه سانتی‌گراد و ۱۶۰۰درجه 
سانتی‌گراد تولید می‌شود. در دماهای بالاتر از ۱۶۰۰درجه 
سانتی گراد واکنش وارد فاز مایع می‌شود و متعاقبا مکانیزم‌های 
متفاوتی درگیر می‌شوند. تشکیل سیلیسیم نیترید در سه مرحله 
سینتیکی رخ می‌دهد. در مرحله اول» نرخ واکنش به‌دلیل 
تشکیل جوانه‌های اولیه سیلیسیم نیترید. تقریبا برابر با صفر 
است. در مرحله دوم واکنش با واکنش سیلیسیم و نیتروژن 
از طریق نفوذ واکنشگر در محل واکنش پیشرفت می‌کند. در 
واقع بیشترین میزان نیتریده‌کردن در این مرحله صورت 
می‌گیرد. در مرحله نهایی. راه‌های نفوذی به داخل هستة 


ذرات با محصولات واکنش بسته می‌شوند؛ بنابراین گاز 


مدل‌سازی اثر مهم‌ترین عوامل موثر بر رفتار سیتتیکی واکنش نیتریده‌کردن... 


واکنشگر باید مسیر جدیدی میان محصولات واکنش بیابد تا 
نیتریده‌کردن قرص کامل شود بنابراین نرخ واکنش کاهش 
می‌یابد. 

در تمام بازة دمایی که واکنش قابل انجام است (دماهای 
پالاکر ات ۱۳۵۸ ره انش راد اکن رتور 
ممانعت ترمودینامیکی برای انجام واکنش وجود ندارد. با 
پیشرفت واکنش. هردو واکنش شیمیایی و نفت‌جوشی در 
پودر سیلیسیم منجر به تغییر در ریزساختار ذرات می‌شود. 
این شرایط چگالی حفرات را تغییر داده و بر سینتیک واکنش 
موثر است. نفوذ موثر گاز واکنشگر مستقیما به حفرات قرص 
وابسته است؛ بنابراین نفوذ موثر گاز واکنشگر با پیشروی 
واکنش تغییر می‌یابد. چندین پارامتر مانند خواص فیزیکی 
قرص سیلیسیم (اندازه. چگالی اندازة حفرات و نوع حفرات) 
و خواص گاز (ترکیب. نرخ سیلان و فشار) می‌تواند روی 
نفوذ واکنشگر مثر باشد. یک مدل موثر که می‌تواند تمام 
پارامترهای اشاره‌شده را فراگیرد. مهم‌ترین مسئله در پیش‌بینی 
رفتار سینتیکی نیتریده‌کردن است. 

تلاش‌های متعددی برای مدل‌سازی واکنش صورت 
گرفته است [7-15]. بیشتر اين مدل‌ها برای توضیح شرایط 
واکنشی گازجامد توسعه داده شده‌اند. [7,16-19]. مدل‌ها در 
اگوی ای تم هر من ره زد من 
ذره‌قرص [19] و مدل واکنش حجمی [16]. مدل سطح 
مشترک تیز (0061 1260عاص1 حتقطه) [7,16] به‌عنوان شر ایط 
اب ام باکت ماهبا رابگر شم ای 
اولیه جامد و غیر متخلخل‌اند. بیان شده است. ایده اصلی پشت 
این فرضیه انجام واکنش در سطح مشترک جامد واکنش‌کرده 
و واکنش‌نکرده است. 

جندین مدل پراساس مدل ٩1۷1‏ توسعه داده شده است. 
چنگ و همکاران [9] مدلی را بر اساس سیستم کنترل نفوذی 
توسعه داده‌اند. آن‌ها فرض کرده‌اند که ذرات داخل قرص 
شرفتشلهل هبل دی مدل پر شعهوآ فک رابظه یس مان 
و نرخ تبدیل به‌صورت تابع خطی بیان شده است که ثابت 


سال سی و دوم شماره دی ۱۶۰۰ 


احسان شاه‌محمدی- فرهاد گلستانی‌فرد- علیرضا می رحبیبی 


همکاران [8] رابطه‌ای را توسعه داده‌اند که مکانیزم‌های موازی 
انجام واکنش و همچنین مکانیزم‌های مختلف انجام واکنش 
در مراحل مختلف سینتیکی را در مدل توسعه‌داده‌شده در نظر 
گرفته‌اند: که و کر کوری: 1201 تأثیر تعیر اندازه دره: در مدل 
ریاضی را در نظر گرفته‌اند. از سوی دیگر حضور هیدروژن 
و اکسیژن در گاز واکنشگر پیش‌بینی سینتیک واکنش را 
پیچیده‌تر می‌کند. بدین منظور درویسبگویچ و رایلی [21] 
مدل دومنطقه‌ای برای واکنش قرص فشرده سیلیسیم در 
اتمسفر نیتروژن /هیدروژن را توسعه داده‌اند [21]. 

عموما مکانیزم‌های متفاوتی در کنترل سینتیک 
نیتریده‌کردن سیلیسیم در حین واکنش حضور دارند. 
مکانیزم‌های درگیر در مراحل اولیه واکنش می‌توانند کاملا 
متفاوت از مکانیزم‌های حاکم بر مرحله نهایی سینتیکی باشند. 
به‌علاوه در زمان‌های به عصوصی از پیشروی واکنش, چندین 
مکانیزم می‌تواند به‌صورت موازی عمل کند. بدین جهت 
مدل‌های اشاره‌شده نمی‌توانند تمام پیچید گی‌های بیان‌شده را 
لحاظ کنند. درنتیجه همیشه اختلاف قابل توجهی بین نتایج 
آزمایشی و پیش‌بینی مدل‌ها در شرایط مختلف واکنشی وجود 
خواهد داشت. از سوی دیگر پارامترهای مژثر زیادی در 
واکتش دخیل هستند که حل دقیق برای همه آن‌ها بسیار 
پیچیده و زمان‌بر خواهد بود. برای حل این موضوع یک روش 
محاسبات نرم بر اساس یادگیری ماشین که بتواند رابطةْ بین 
ورودی‌ها و هدف مسئله را بیابده کاربردی بوده و چنین 
روشی دارای مزایای بسیار زیادی در مواجهه با چنین شرایط 
پیچیده و غیرخطی است. 

در دهه‌های اخیر. نظریات یادگیری ماشین مانند شبکة 
عصبی مصنوعی 2600216 لحسعه لممتیه) للم و 
5 به‌عنوان رایج‌ترین روش‌های محاسبات نرم برای 
حل‌کردن چنین محدودیت‌ها و پیچیدگی‌های محاسباتی به 
کار گرفته شده‌اند [22-25]. نتایج نشان‌دهندة عملکرد موفق 
نظریات داده کاوشی است. به‌هرحال. مدل‌های ۸ و 
5 دیدی از مدل ایجادشده نمی‌دهند و به‌راحتی 


فرمول‌های تجربی نمی‌توان از آن‌ها استفاده کرد. در میان 


سال سی و دوم شماره دو ۳۹ 
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روش محاسباتی نرم. شبکة 01۷ به‌عنوان روش 
خودسازمان‌دهنده (81]0782201260) برای مدل‌سازی و 
کشف رفتار سیستم‌های پیچیده و ناشناخته با استفاده از 
داده‌های ورودی و خروجی شناخته شده است (ایواحننکو 
[26,27]). هدف اصلی این مطالعه بررسی کارآمدی شبکه 
11 برای پیش‌بینی سینتیک نبتریده‌کردن سیلیسیم است. 
مزیت اصلی روش ]0۷ در مقایسه با سایر روش‌های 
مشابه مثل ۸ تعیین روابط بین داده‌های ورودی و 
خروجی به‌صورت پارامتری در قالب یک رابطه بیان می‌شود. 
درحالی‌که در روش ۸ این روابط به‌صورت پنهان 
درون‌ساختار شبکة مصنوعی است. برای توسعة مدل کارآمد 
بر اساس 01۷19۳1 داده‌های جامعی از مقالات شامل ۲۱۸۲ 
نتیجهٌ آزمایشی به کار گرفته شد. مدل 01/۳۳ 
توسعه‌داده‌شده ارتباط بین درصد تبدیل سیلیسیم به نیترید را 
با زمان دما؛ درصد نیتروژن گاز واکنشگر انداز قرص و 
اندازءٌ ذره سیلیسیم در قالب روابط بیان می‌کند. نتایج 
به‌دست‌آمده از 20111) با دو مدل عمومی موجود دراین‌زمینه 
از طریق شاخحص‌های خطای آماری مقایسه می‌شود. اهمیت 
نسبی پارامترهای موثر بر درصد تبدیل سیلیسیم به نیترید با 
استفاده از آنالیز حساسیت بررسی می‌شود. قدرت مدل 


1۳11 ارائه‌شده با تحلیل پارامتری تأیید می‌شود. 


روش تحقیق 
11 به‌عنوان یک روش یادگیری ماشین برای تشخیص 
الگوهای غیرخطی پنهان‌شده درون داده‌ها معرفی شد. این 
روش اولین‌بار از سوی ایواخننکو در ۱۹۷۱ معرفی شد [26]. 
رابطةٌ غیرخطی بین متغیرهای ورودی و جواب‌ها به‌صورت 
توابع چندجمله‌ای ارائه شد. در حقیقت. الگوریتم بهترین 
چندجمله‌ای درجه دوم را با تعیین بهترین ضرایب برای 
لحاظ کردن چندجمله‌ای به‌عنوان وابستگی متغیرها در هرلایه 
انتخاب می‌کند. الگوریتم 01۷1011 دو دسته داد آموزشی و 
صحت‌سنجی را برای توسعة مدل به کار برده است و بهترین 


خروجی‌ها براساس بهترین ترکیب از داده‌های ورودی مستفل 


پیچیده‌ای از 


به‌دست می‌آید. برای دستیابی به اين ام ن 
سری توابع ولترا (۷۵1]۵02) به‌عنوان چندجمله‌ای ولترا- 
کلمگرو به‌کاررفته است تا بهترین ارتباط بین مقادیر پاسخ و 
متغیرهای چندورودی به‌صورت رابطهٌ (۱) تعیین کند. 
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به‌طوری که ,۱ تعداد متغیرهای م ستقل ا ست. ماتریس 
داده‌های ورودی که شامل ,لا متغیر مستفل و لا جواب 
مشاهده شده است. در رابطه (۲) ارائه شده است. ماتریس 


سمت چپ شامل جواب‌های مشاهده شده و سمت راست 


شامل ,لا متغیر مستقل است. 
سس حول وی 3 
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جندجمله‌ای درجهة دوم برای متغیر مسستقل در رابطةٌ 
(۳) ارائه شده آمتتتتا: برای بل[ متغیر رادطة در حه دو 


2 


2 
و25 مره 1 


وز*م 23 ‌ِ 220 ‌ِ 210 ِ 0 < ۷ 


۳( 
با استفاده از روش رگرسیون ضرایب چندجمله‌ای 
درجه دوم تعیین می‌شود. برای دستیابی به بهترین تطابق بین 
رش ال و داده ها واقی: ووشن اقا فان مررعازت 
| ستفاده می شود. هدف | صلی روش 0۷177 حداقل کردن 
مربع تفاضل بین داده های خروجی و مقادیر واقعی طبق 

رابطةه (۶) است. 


2 فا 
( ۶:۱ ]2 < ‌ 
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به‌طوری که داده پیش‌بینی شده است. 


مدل‌سازی اثر مهم‌ترین عوامل موثر بر رفتار سیتتیکی واکنش نیتریده‌کردن... 


در این مطالعه یک مدل پیش بینی کذنده بر اساس روش 
1 توسعه داده شده است تا درصد تبدیل سیلیسیم به 
نیترید را در رابطهٌ (۵) تخمین بزند. 


(6( رم لا ولد ی + و1٩3‏ 


8) 2 
برای توسعه مدل» یک دسته داده جامع برگرفته از نتایج 
آزمایش‌های چنگ و همکاران استفاده شد. مجموع داده‌ها 
شامل ۲۱۸۰ نقاط داده است. چندین پارامتر موثر شامل زمان 
دماء اندازهُ قرص.ء درصد پیتروژن و اندازه ذرات سیلیسیم 
به‌عنوان داده‌های ورودی در نظر گرفته می‌شوند. اين پارامترها 
به‌عنوان مهم‌ترین پارامترهای موثر در نفوذپذیری گاز 
واکنشگر شناخته می‌شوند. واکنش نیتریده‌کردن سیلیسیم از 
نظر سینتیکی در دسته واکنش‌های کنترل نفوذی طبقه‌بندی 
می‌شوند که سینتیک آن‌ها به نفوذ گاز واکنشگر وابسته است؛ 
بنابراین» پارامترهای نشان‌داده‌شده به‌عنوان متغیرهای 
پیش‌بینی کننده برای تخمین میزان تبدیل به نیترید انتخاب 
می‌شوند. بازهُ متغیرهای ورودی و خروجی در جدول (۱) 
ارائه شده است. شاخص‌های آماری شامل حداقل, حداکش 
میانگین و انحراف معیار نیز برای هر پارامتر در جدول (۱) 
آمده است. باز؛ُ انتخابی برای متغیرهای ورودی. بازه‌های 
عملی (واکنش اغلب در این بازه‌ها رخ می‌دهد) پوشش 
می‌دهد. 
برای توضیح بیشتر, نمودار میله‌ای متغیرهای ورودی و 
خروجی برای بازه‌های متفاوت در شکل (۱) نشان داده شده 
است. باید متذکر شد که مدل ارائه‌شده در بازه‌هایی که داده‌ها 
تمرکز بیشتری دارند قابل اعتمادتر است. در این مطالعه؛ یک 
پارامتر جدید معرفی شده است تا اثر اندازه ذرٌ سیلیسیم را 
لحاظ کند. این پارامتر شامل ۵ بازه از اندازه ذرات به‌صورت 


حدول (۲) است. 


سال سی و دوم شماره دو ۱۶۰۰ 


احسان شاه‌محمدی- فرهاد گلستانی‌فرد- علیرضا میرحبیبی ۳۱ 


جدول ۱ متغیرهای موثر و داده‌های آماری 

متغیر نماد | حداقل | حداکثر | میانگین | انحراف معیار 
زمان [صنص] گ | هی | ۲۳۹۸ ] ۱۰۱,۸۷ ۷/۳۸۳ 
دما [>1] ۴ | ۱۵۵۸ | ۱2۵۸ | ۱۰۶,۱۲ ٩1۸2‏ 
انداز؛ قرص [7070] ۳ ۲ ,۵ ۲,۵ 
درصد نیتروژن [96] ,۱ ۱.۰ ۸۵,۷ ۲۲,۱ 
اندازة ذرات سیلیسیم[] | ٩‏ ۱ ۵ ۱۶ ۱۵ 
درصد تبدیل به نیتروژن [96] | * | ۰,۱۱۷ | ۳,۶ | ۱۵,۹۷ ۹۸ 


1800 
300 
1600 
1400 250 
1200 
5 2۳200 2 
0 شِ 5 5 
و و 130 3 3 
۳ 10 ۳ 
40 
50 
200 
0 0 
8 1618 1598 1578 1558 1538 1518 1498 1478 1458 40 30 20 10 0 250 200 150 100 50 0 
[6] ستعهومه1 [۲۵] ٩‏ ۵۶ جمزجه:ه60 [منس] ] 
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1500 1500 1500 

<< حر حر 

4 4 4 

ب ۴ 

و و دِ 
1000 2 1000 2 1000 2 
500 500 500 
0 0 0 

5 4 3 2 1 100 80 60 40 20 0 12 10 8 6 4 2 
6 16ءنانعظ جم‌ناز۹ [۷۵] ۳6۱6۵۵1286 طعع00ز( [حصصط] عم2رو ۳۵۱۱6۲ 
شکل ۱ نمودار متغیرهای ورودی و خروجی 
حدول ۲ توزیع اندازهٌ ذرات 
5و۳ ۳9-4 3و۳ ۳922 ۳-1 


> ص63 | هب ۳5>63> صن 53 | هی ۳553 > یرک[ صیر ک۳5>4> هن 38ص 23>38 


برهم‌نهی داده یکی از مشکلات مهمی است که در حین تقسیم شدند. داده‌های آموزشی برای توسعدٌ مدل جدید 
آموزش الگو با استفاده از روش یادگیری ماشین رخ می‌دهد. پایه‌ای استفاده می‌شود و داده‌های آزمایشی برای آزمودن 
بنژاف و همکاران در سال ۱۹۹۸ روش جدیدی را برای عملکرد مدل توسعه‌یافته به کار می‌رود. برای دستیابی به 
کاستن این برهم‌نهی و بهبود عمومیتبخشی مدل توسعه‌یافته  .‏ بهترین و عمومی‌ترین مدل. بایستی عملکرد مدل آموزش یافته 
معرفی کردند. در این روش داده‌های ورودی اولیه به‌صورت . برای هردو دسته داده آموزشی و آزمایشی یکسان باشد. برای 
اتفاقی به دو دسته داده‌های آموزشی و داده‌های آزمایش دستیابی به مدل قابل اعتماه حداکش حداقل. متوسط و 


سال سی و دوم شماره دی ۱۶۰۰ نشریه مهندسی متالورژی و مواد 


۳۲ 


انحراف معیار پارامترها بایستی به پایداری مقداری در 
داده‌های آموزشی و آزمایشی برسند. در این مطالعه ۲۱۸۰ 
نقطه داده برای توسعهة مدل نهایی استفاده شده است؛ بنابراین» 
۰ درصد داده‌ها (۱۷۶۹ نقطه داده) برای آموزش به کار رفته 
و ۲۰ درصد باقی‌مانده (۶۳۷ نقطه داده) برای آزمایش به کار 
رفته‌اند. 

بعد از تقسیم‌بندی داده‌هاء داده‌های آموزشی به الگوریتم 
211 اعمال شده و روابط به‌دست‌آمده از مرحلة آموزش 
مدل برای تخمین درصد پیشرفت نیتریده‌کردن در مجموعةٌ 


روابط )1 ال اشتتا: 


4(۰حوممتنهآظ ۵۶ مامت 1 12۵۲ 


- 10751۲ *< 3.7 - 0.25 + ]0.29 + 260- < م6 
2- 107512« 8.1 


)2 
۲۳ + 0.111 + 0.057۳ + 2.9- < رز 
2 - 0.0002(]۳2 + 
(10) 
- 0.00385 - ]0.19 + 42۳۴5 + 0.47 < و1 
12 - 0.81۳52 
)0 
+ ]0.00073 + 0.03271 + ۱۳ - 3.7 < ور 
72 - 107511۳2 *< 9.1 
)0 
۰(عفصمتنهآظ 0۶ مامت ۳ 2 12۷۵۲ 
0.08126 + 0.7266 + 0.552 + 1.7- < ورزر 
22 - 0.04323 
13 


+ و0064 + مل(0.33 + و0516 - 1.4- < ور 
08 -- 0.0261 


9 
ولاع0.06926 + 0.42266 + 0.5269 + 4.2- < ور 
12 - 0.02563 
(2( 
:(2عفص ۱۵0 0۶ عمجت ) 3 12۵۲ 
+ و031 + ور(0.89 + 0.83- < ۸.6 
وگل0.0067 - 0.061 - ورلار 0.061 
9 


مدل‌سازی اثر مهم‌ترین عوامل موثر بر رفتار سینتیکی واکنش نیتریلهکردن... 


ِ ر6 0.0542 ٍِِ۳ 9111 02 ِ و072 ِ 1 ت 7 
0.0088 + 0.038 


(1) 
۰( 1حوممتنهآظ ۵۶ مامت ۵۲ 4 وه 
و14 + م0035 - 0.61- < (9۵) ۳ 
م12 ۲ ر0.382 + میالم 0.792 
۹1 


شکل ۲ سیلیسیم خالص فشرده پس از پرس پودر 


روش آزمایش 
مادهُْ به‌کاررفته در این تحقیق پودر سیلیسیم با خلوص 
0٩‏ درصد ( فصن ...رامعم ...هو -عمجر1[ 


«س476-و)_ است که به‌صورت کمی قبل از انجام 
آزمایش‌ها مشخصه‌یابی شده است. سایر عناصر اضافی پودر 
شامل 0۵ /۰/۰۰۰۱> وع11 /۰/۰۰۰۲> است. به‌علاوه 
انداز؛ ذرات پودر سیلیسیم با استفاده از آنالیز انداز دانه 
(1۷۲۸۲,۷۲ ,۷1۸1,100229) تعیین شد. پودر اولیه در قالب 
فلزی با فشاره۷1۴ ۰,۱ فشرده و فرص سیلیسیمی ایجاد شد 
(شکل (۲)). قرص فشرده سپس به دانه‌های کوچک‌تر 
شکسته شد تا فضای تماس با بوته به حداقل برسد. 

باتوجه‌به نتایج مدل برای زمان اتمام واکنش. ۰آبازه 
زمانی برای محاسبهٌ درصد نیتریده‌کردن در ۱۵۷۳ کلوین 
انتخاب شد. 


به منظور انجام فرایند نیتریده‌کردن سیلیسیم کوره تیوبی 


سال سی و دوم شماره دی ۱۶۰۰ 


احسان شاه‌محمدی- فرهاد گلستانی‌فرد- علیرضا می رحبیبی 


از جنس آلومینا استفاده شد. از آنجایی که اکسیژن با ایجاد 
لاه اکسیدی نازک در سطح سیلیسیم. واکنش نیتریده‌کردن را 
به تعویق می‌اندازد. مقدار کمی گاز هیدروژن به گاز نیتروژن 
اشتافه کل 9ا ‏ زد کی افبارمهته ۵ تیار اکشفن اخیاان 
موجود در گاز واکنشگر را حذف کند. بدین جهت ترکیب 
حجمی ۵درصد هیدروژن-۹۵درصد نیتروژن به‌عنوان گاز 
واکنشگر در حین فرایند استفاده شد. نمونه‌ها در کوره قرار 
گرفتند و به‌مدت تعیین‌شده در معرض گاز واکنشگر در دمای 
که وهی ی ای ۱ باه ای 
متفاوت انجام گرفت و میزان پیشروی واکنش با افزايش وزن 


نمونه‌ها اندازه‌گیری شد. 


انداز؛ درصد ثیتریده کردن 
به‌دلیل اينکه وزن ذرات سیلیسیم بعد از واکنش با گاز نیتروژن 
و تشکیل نیترید سیلیسیم افزايش می‌یابد. روابط 
استوکیومتریک بر اساس رابطة (۷) وابستگی بین درصد تبدیل 
نیتروژن را با افزايش وزن بیان می‌کند. همچنین نتایج 
آزمایشگاهی در جدول (۳) خلاصه شده است. 


2 ملاین < ول( 2+ 35 
)۸( (م)ملاونگ ما رل( ؟ه )دوه ۷۷ < ۸۷۷ 
)0( «ورعامط ع رام < ک 
۰ مامح < 3. ۳ 
۱۱۱ روم - 9 (3. ۳ 


در روابط نشان‌داده‌شده )۷ تا ۲ 2 بیانگر درصد 
واکنش» ۸ افزایش وزن. ۷۷0 وزن اولية قرص. زع1۷ وزن 


مولکولی سیلیسیم و 2« وزن مولکولی نیتروژن است. 


سال سی و دوم شماره دو ۱۶۰۰ 


۳۳ 


۹ ۸۳۷ ۷۳۷ زمان (دقیقه) 
۱۱۳۰ ۳۳ 1۳۸ ۳ 
۱۹:۹ ۰۳۱ .پآ ۲۵ 
۱9,۶ ۰۳۱ ۷«رپآ« ۰ ۳۷ 
۱9۹,۸3 ۰۳۲ ۸آ ۵۰ 
۱۸,۶7 11 ۰۳۲۵ ۹ 
۱۸,۹ ۵ اِ/۳ ۷۵ 
۱۸۹۷/۵ ۳۸« ءٍرُآ: ۸۷ 
۳۰۰1 11 ۳۹۹ ۱ 
۱۹,۹۲ ۳« ۲۷۱ ۱ 
۳۱,۲۳ ۸ ۰۳۳۹ ۱۲۵ 
۲۰۶۱ ۰9۲ بآ ۱۳۷ 
۱۹/۳ ۲" 15۲ ۱۰ 
۳۲۳۸ ۲ ۱۳۶ ۱3 
۳۱۳,9۵ 11 ۳۹3 ۱۷۵ 
۲۱۷ ۵ ۳۱ ۱۸۷ 
۲۳,2۵ ۰ ۷"رصپآ« ۰ ۲ 
۲۳,۹2 9« ۳۳۸ ۳۲ 
۳۳/۸۳ ۹ ۰,۲۸۹ ۳۳۵ 
۲۳,۸ ۵" ۵أآ ۰ ۳۳۷ 
۱:9۳ ۲ ۰۳۷۹ ۲۵ 

نتایج و بحث 


در این مطالعه. سینتیک تشکیل بلاوله با استفاده از 
آزمایش‌های تجربی و نظریةٌ محاسبات نرم بررسی و نتایج 
برای صحت‌سنجی مدل پیشنهادی با یکدیگر مقایسه و تطابق 
قابل قبول مشاهده شد. به‌علاوه مدل صحت‌سنجی‌شده با 
مشاهدات آزمایشی سایر مدل‌های به‌دست‌آمده در تحقیقات 
مشابه مفایسه شد. 

نتایج الگوی اشعةٌ ایکس (680 و تصویر 
میکروسکوپ الکترونی روبشی. پیشروی واکنش و 
مورفولوژی فازهای تشکیل‌شونده را در گام‌های متوالی با 
گذر زمان در دمای ثابت نشان می‌دهند. علاوه بر اين» عملکرد 


مدل توسعه‌یافته با پارامترهای آماری ارزیابی شد. سپس 


۳ 


مهم‌ترین پارامترهای تخمین درصد تبدیل با استفاده از آنالیز 
حساسیت تعیین می‌شود. برای صحت‌سنجی بیشتر مدل, 
آنالیز پارامتریک صورت گرفت. این آنالیز از تبعیت نتایج مدل 
با واقعیت فیزیکی حاکم در مسئله اطمینان حاصل می‌کند. 


آنالیز ریزساختار 
ریزساختار قرص نیتریده در ۱۵۷۳ کلوین بعد از آماده‌سازی 
اولیه با استفاده از 875186 مشخصه‌یابی شد (شکل (۳)). تشکیل 
فاز ۵51314 طی مکانیزم بخار مایم جامد (16) یا 
به‌صورت ساده با فرایند رسوب شیمیایی بخار (۷) بین 
سیلیسیم گازی و نیتروژن رخ می‌دهد. بیشترین درصد تشکیل 


۷۶6۸۱ ۸۵۷ 
‌ 


1۲۵0 ۱0۲۱, 861. ۰ 


مدل‌سازی اثر مهم‌ترین عوامل موثر بر رفتار سینتیکی واکنش نیتریلهکردن... 


واکنش گازجامد در سطح سیلیسیم به‌دلیل نفوذ اتمی نیتروژن 
درون فاز یا در حضور فاز مایع که تشکیل این فاز را تشویق 
کناتساد تراد 

الگوی پراش اشعهٌ ایکس محصولات آشکار می‌کند که 
قرص سیلیسیم واکنش کرده شامل سیلیسیم واکنش‌نکرده. 0 
ولگ و بل8-51 بوده و از سوی دیگر پیک اکسید سیلیسیم 
(«5(0) و سیلیسیم اکسی‌نیترید (5(2:0) مشاهده نشدند. 
تسبیت. فاز 0 مشاهد‌شده سیان بشتر از فاز.9] است: 
همان‌گونه که در شکل (1) دیده می‌شود با پیشروی واکنش 
شدت پیک سیلیسیم کاسته می‌شود و از سوی دیگر شدت 


پیک هردو فاز 0 و افزايش می‌یابد. 


۱۷۸/۵۰ 10.26 ۲ 
06۱: 6 
0316)۳0/0۸(: 9 


8۶0 ۳۱۷: 0,0 ۷ 
۱۷۱۵۷۷ ۵۱۵: 21.67 ۷ 
56۱ ۷۸: 0,0 


شکل ۳ ریزساختار قرص نیتریده در دمای ۱۵۷۳ درجه کلوین 


نشریه مهندسی متالورژی و مواد 


سال سی و دوم. شماره دی ۱۶۰۰ 


احسان شاه‌محمدی- فرهاد گلستانی_فرد- علیرضا می رحبیبی 
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شکل ۶ الگوی 0 مسه قرص که به‌ترتیب در سه زمان ۲۵۰6۵ ۱۰۰0 6) ۲۵ دقيقه در دمای ۱۵۷۳ درجه کلوین نیتریده شده‌اند 


آنالیز عملکرد 
برای تعیین عملکرد مدل توسعه يافتةٌ 113 ۶ پارامتر 
خحطای آماری شامل بایاس؛ ریشه دوم میانگین خطا 


رفته است. 


(00-:۳) ی 
آ 


,۱۳ < 125ظ 


3 02 - کر - 8۱۸58 


(م0-:0)(مو۳-ز۳) یز - 
:۱/60( حرط زظ) لو 


(۱۵( 


۳2 -0) ی 


بت 1 << 12 
0۱-0۶ 


۱۹0( 
به‌طوری که نل) م قدار مسا هده‌شده؛ :ظ م قدار 

پیش بینی شده به‌وسیلة مدل» تعداد نقاط داده, 0 مقدار 
متو سط مشاهدات و 3 مقدار متو سط مقادیر پیش‌بینی شده 


است. 


سال سی و دوم شماره دو ۱۶۰۰ 


نتایج مدل 0۷]«۳) توسعه‌داده‌شده و مشاهدات 


آزمایشگاهی واقعی برای آموزش و آزمایش در شکل (۵) 


() و (0) نشان داده شده و به‌علاوه خطاها و بهترین توزیع 
نرمال روی این خطاها در این شسکل نشان داده شده‌اند. 
تطابق خوبی بین درصد تبدیل پیش‌بینی‌شده با مدل 01۷1۲۳1 
و داده های واقعی وجود دارد. اين امر بیانگر این واقعیت 
است که عمومی سازی مدل تو سحه‌داده شده برای داده‌های 
تصادفی در بازهُ آموزشی مدل قابل قبول است. از سوی 
دیگر. خطاهای یک مدل پیش‌بینی‌کننده خوب بای ستی روند 
و توزیع خاصی را دذ بال 5ند [28]. اين ویژگی برای 
لحاظ کردن عدم‌قطعیت در طراحی مدل پیش‌بینی‌کننده قابل 
استفاده است. همان‌گونه که در شکل (۵- 2 و 0) نشان داده 
شده است. خطاهای مدل 01۳1 توسعه‌یافته تقریبا توزیع 
نرمالی را برای هردو د سته داده آموز شی و آزمای شی دنبال 
می‌گفد؛ ب نابراین. عدم قطء‌ یت مدل 01۳01 می توا ند 
به‌سادگی با استفاده از لحاظ کردن پارامتر اتفاقی با یک توزیع 


نرمال محاسبه شود. 


۳۹ 


مدل‌سازی اثر مهم‌ترین عوامل موثر بر رفتار سینتیکی واکنش نیتریده‌کردن... 
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شکل ۵ مقایسه نتایج پیش‌بینی مدل و مشاهدات آزمایشگاهی ه)آموزشی 0) آزمایشی 


برای محاسبه کمی عملکرد مدل توس عه یافته 
پارامترهای خطای آماری در جدول (۶) برای هر دو دسسته 
داده آموزشی و آزمایشی ارائه شده است. درواقع. پارامتر *1 
بیانگر ضریب همبستگی بین دو مقدار اندازه‌گیری‌شده و 
بش اه اسستیتو ا گر باواهش. بت اه ماش هه 


تشبسان داهتلم: تیک قوعع بیرخ مقادان پیت مش تسده 


شرب مچنلسی متالزرژی و مراد 


شاه شاه اسفی هی اه عطالهه یی یی نها 
شده است. زیرا ارزیابی صریح‌تری از عملکرد مدل را بیان 
می‌کند. به‌علاوه پارامترهای بایاس و 1815۳ در جواب‌های 
دقیق‌تر مقادیر نزدیک به صفر دارند. مطابق جدول (4) مدل 


171 عملکرد یکسانی برای هر دو دسته داد آموزشی و 


آزماٍ شی دارد. مقادیر بایاس و 80۷1517 یک سان و نزدیک به 


سال سی و دوم شماره دی ۳۹2 


احسان شاه‌محمدی- فرهاد گلستانی_فرد- علیرضا می رحبیبی 


صفر هستند. مدل 0۳1 با مقادیر خطای ۰/۹11- و 
۵ و مقدار 12۶ غریبا برابر با ۰/۹۳ برای داده ها 
آ ویو میتی معارهای وتات رب 
همبستگی قوی را تحقق می‌بخشد. برای توضیح بیشتر 
همیستگی بین خروجی‌های مشاهده شده و پیش‌بینی شده با 
مدل, مقایسه‌ای بین داده‌ها هم برای داده‌های آموزشی و هم 


۳۷ 


آزمایشی به‌همراه میزان خطا در شکل(1) آمده است. در این 
شکل بهترین رگرسیون خطی بین داده‌های واقعی و خروجی 
مدل نشان داده شده و معادلهة حطی برازش نزدیک به 
معادلة بهينة ۷-۶ است. نتایج بیانگر این حقیقت است که 
مدل توسعهیافته توانایی پیش‌بینی مناسب و قابل‌اعتمادی 


دارد. 


جدول ۶ صحت‌سنجی مدل پیشنهادی با استفاده از معیارهای خطای آماری 


داده‌های آزمایشی داده‌های آموزشی شرایط رابطه نوع پارامتر خطا 
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شکل 1 مقايسة بین خروجی پیش‌بینی‌شده و نتایج آزمایشگاهی 


سال سی و دوم شماره دو ۱۶۰۰ 


۳۸ 


برای ارزیابی بیشتر عملکرد مدل 01011 توسعه‌یافته. 
پیش‌بینی این مدل با دو مدل قبلی لی و همکاران (۱۹۹۷) و 
چنگ و همکاران (۲۰۰۰) مقایسه شده است. مدل لی و چنگ 
دو مدل شناخته‌شده برای پیش‌بینی سینتیک نیتریده‌کردن 
سیلیسیم به‌روش نیتریده‌کردن مستقیم است و چندین پارامتر 
موثر و مکانیزم که پیروی واکنش را کنترل می‌کند در آن‌ها 
لحاظ شده است؛ بنابراین مدل‌های اشاره‌شده می‌توانند 
به‌عنوان مرجع قابل قبول برای ارزیابی مدل 01011 
توسعه‌یافته به کار روند. برای مقايسة عملکرد مدل دما و 
درصد گاز نیتروژن در اتمسفر به‌صورت ثابت به‌ترتیب ۱۵۷۳ 
کلوین و ۹۵درصد در نظر گرفته شده‌اند. این مقادیر برای دما 
و غلظت گاز. عمومی‌ترین حالت درنظرگرفته‌شده برای انجام 
واکنش است که اغلب برای اطمینان از انجام واکنش اعمال 
می‌شود. براساس مقایسه‌هایی که بین نتایج دو مدل نشان داده 
شده است (مدل 01۷9۳1 توسعه‌یافته و مقادیر آزمایشی 
داقتی 1 تمودان انت‌ستتی. فزسیم. شد. هکل ۰ (): 
همان‌طور که نشان داده شده است. هر دو مدل لی و چنگ 
در ابتدای واکنش پیش‌بینی‌های خوبی داشتند و نتایج با 
داده‌های آزمایشی واقعی تطابق دارند ولی بعد از دقیقه ۱۷۵م 
پیش‌بینی این مدل‌ها دارای خطای قابل توجهی می‌شود. 
مطابق این شکل مدل توسعه‌یافته تطابق قابل ملاحظه‌ای با 
داده‌های آزمایشگاهی در تمام مراحل انجام واکنش را داشته 
و ضریب خطای 132-۰0/۹۳ این امر را تأیید می‌کند. درحالی که 
مدل لی با خطای ۰/۱۳ و چنگ با ۰/۵7 پیش‌بینی خوبی 


نداشتند. 


اس 
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شکل ۷ مقایسه بین نتایج مدل لی(۱۹۹۷) و چنگ(۲۰۰۰) و مدل 
11 توسعه‌یافته با مقادیر آزمایشگاهی واقعی 


مدل‌سازی اثر مهم‌ترین عوامل موثر بر رفتار سینتیکی واکنش نیتریده‌کردن... 


آنالیز حساسیت 

برای تعیین مهم‌ترین پارامتر در پیش‌بینی درصد پیشروی 
واکنش آنالیز حساسیت براساس آزمون گاما انجام شد. آزمون 
گاما روشی برای تعیین درجه اثر ورودی روی خروجی مدل 
تاه ارت آزمون در مدل‌هایی با چندین ورودی صورت 
می‌گیرد. انتخاب ترکیبی از مهم‌ترین متغیرها یکی از مسائل 
بسیار مهم در مدل‌سازی فرایند با روش‌های یادگیری ماشین 
مثل 7 است [30]. 


برای تعیین میزان اهمیت و اثر پارامترهای ورودی» شش 
سناریو لحاظ شده است. در سناریو اول. تمام پارامترهای 
ورودی در آنالیز گاما لحاظ شده‌اند. در سناریو بعدی» تمام 
پارامترهای ورودی تکبه‌تک از آنالیز گاما خارج می‌شوند. 
تغییرات در پارامترهای گاما مثل وزیمم۷ بعد از حذف پارامتر 
محاسبه می‌شود. تغیبرات بیشتر در وززع۷ نشان‌دهنده اهمیت 
بالای پارامتر حذف‌شده در تخمین درصد تبدیل به نیترید 
انسگ: وبیر۷ تالا تشان‌دهنده اعمی الا باران عذف شاه 
بوده و انحراف بیشتری در پیش‌بینی‌ها نسبت به مقادیر واقعی 
ایجاد می‌کند. 

نتایج آنالیز گاما در جدول (۵) خلاصه شده است. برای 
توضیح بیشتر نتایج آنالیز گاماء مقادیر نیعم برای سناریوهای 
مختلف در شکل (۸) نشان داده شده است. طبق جدول (۵) 
زمان. درصد نیتروژن گاز واکنشگر و دما به‌ترتیب 
اثر گذارترین پارامترها در تخمین درصد تبدیل نیتروژن است. 
در حقیقت این نتایج با واقعیت‌های فیزیکی مسئله تطابق 
دارد. عموماء سیلیسیم نیترید از دو فاز بلالوگ-ه و بلامنک-8 
تشکیل‌شده است [31]. تشکیل دو فاز در دو مکانیزم نفوذی 
مجزا رخ می‌دهد که اين مکانیزم‌ها کند بوده و تحت‌تآثیر 
متغیرهای متفاوت هستئدة بتایراین: زمان واکنش نقش مهمی 
فد میک فت راکش فارم. که به وک رها این شاشی ۱۶ 
می‌دهد تا درون لایه محصول اطراف ذره نفوذ کنند؛ بنابراین 
آنالیز گاما به‌طور صحیح زمان را مهم‌ترین عامل موثر بر 
پیشروی واکنش تشخیص داده است. 


سال سی و دوم شماره دی ۱۶۰۰ 


احسان شاه‌محمدی- فرهاد گلستانی_فرد- علیرضا می رحبیبی 


جدول ۵ نتایج آنالیز گاما 
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شکل ۸ نتایج آنالیز گاما (مقدار پارامتر ونم برای سناریوهای متفاوت) 


آنالیز 071 همچنین نشان می‌دهد که درصد نیتروژن 
دومین پارامتر موثر روی سینتیک نیتریده‌کردن است. درواقع؛ 
چگالی جوانه‌های سیلیسیم نیترید روی سطح مشترک وابسته 
به میزان نیتروژن جذب‌شده در سطح است. میزان نیتروژن 
جذب‌شده متأثر از فشار نیتروژن روی سطح مشترک واکنش 
است. با افزايش فشان میزان نیتروژن جذب‌شده در سطح 
مشترک افزایش می‌یابد. فشار بالای نیتروژن جوانه‌های اولیه 
اوه بیشتری تشکیل می‌دهد و سینتیک اولیه را افزایش 
می‌دهد. به‌علاوم. فشار بالا مکانیزم‌های نفوذ را با افزایش 
غلظت نیتروژن در حفرات درون ذره‌ای سریع‌تر می‌کند. 

نفوذ که نقش کنترل‌کننده روی نرخ واکنش نیتریده‌کردن 


دارفه شین شمان سار ان ففاش فیح تفر فاگ 


سال سی و دوم شماره دی ۱۶۰۰ 


به‌دلیل انرژی فعال‌سازی بالای آن به‌شدت به دمای متغیر 
وابسته است. مطالعات قبلی [8-10,14,32,33] نشان می‌دهد 
در صورتی نفوذ نیتروژن درون لایه محصول امکان‌پذیر است 
که نیتروژن به‌حالت اتمی فعال باشد. اگر اکسپیژنی نباشد. 
سطوح داغ تنها منابعی هستند که نیتروژن اتمی تولید می‌کنند. 
به‌دلیل اينکه تجزیه نیتروژن به‌حالت اتمی مقادیر زیادی از 
انرژی را می‌طلبد. نفوذ شدیدا به دما وابسته است. از سوی 
دیگر: این امر ممکن است واکنش در حالت بخار صورت 
گیرد. در این حالت بخار سیلیسیم با نیتروژن واکنش می‌دهد 
و دما نقش اصلی را در تولید فشار بخار کافی سیلیسیم ایفا 
می‌کند. نتایج آنالیز گاما در راستای نتایج ثبت‌شده در مقالات 


قبلی و آزمایشات 7 تعیین نرخ واکنش در تحقیق حاضر است. 


آنالیز پارامتریک 
پارامتریک انجام شد تا اينکه از هم‌راستایی مدل 0111 با 
مفاهیم فیزیکی سینتیک واکنش اطمینان حاصل آید. روش 
آنالیز پارامتریک بدین شیوه است که با استفاده از تغیبر یکی 
از پارامترها درحالی که سایر پارامترها در محدوده میانگین 
خود ثابت نگه‌داشته شده‌اند» تغییرات جواب‌های مدل ثبت و 
بررسی می‌شوند. همان‌طور که در بخش قبلی بیان شد. زمان؛ 
درصد نیتروژن و دما از مهم‌ترین پارامترهای مثر بر سینتیک 
واکنش هستند؛ بنابراین تأثیر این پارامترها با به‌کارگیری آنالیز 
پارامتریک بیشتر بررسی شد. نتایج آنالیز پارامتریک 
در شکل 69 ارائه شده 
نیتریده‌کردن با زمان واکنش برای ؛ دمای مختلف 
(۲-1480,1510,1550,1600016) و ۶ غلظت نیتروژن 
متفاوت (10066 ,9560 ,9096 ,2-8590) به‌ترتیب در شکل 
0(۹) نشان داده شده است. باتوجه‌به شکل (۵ با 


است. تغییرات درصد 


افزایش دما و غلظت نیتروژن. درصد تبدیل به نیترید افزايش 
می‌یابد. به‌ازای هر درجه تغییر دمای واکنش ۰/۰۵۵ درصد و 
به‌ازای هر ۱درصد افزايش غلظت گاز نیتروژن ۰/۶درصد 
افزايش پیشروی واکنش در زمان یکسان رخ می‌دهد. واکنش 
نیتریده‌کردن سبلیسیم یک واکنش کنترل نفوذی است که 
به‌شدت به دمای واکنش و غلظت واکنشگرها وابسته است. 
بنابراین با افزايش دما و درصد نیتروژن. نرخ واکنش افزایش 
می‌یابد. به‌علاوه در مراحل اولیةٌ واکنش. سیلیسیم به‌طور 
مستقیم در معرض گاز واکنشگر قرار می‌گیرد درحالی‌که با 
ادامه واکنش لایة محصولات سطح نمونه را پوشانده و پدیده 
نفوذ واکنشگرها مهم‌تر می‌شود. در نتیجه. دما و پارامتر 
غلظت نیتروژن در مراحل بعدی واکنش اثر مهم‌تری دارند. 
باتوجه‌به شکل, نتایج نشان می‌دهند که مدل 6/0 
توسعه‌یافته با الگوهای فیزیکی تطابق دارد. بایستی یادآوری 
شود که تغییرات غلظت نیتریده‌کردن اثر بیشتری در مقایسه 
با دما روی تغییر درصد تبدیل سیلیسیم دارد که این واقعیت 
هم در آنالیز حساسیت و هم در کارهای گذشته بیان شده 
است [8-10] و [13-15] و [20,32,33]. 


مدل‌سازی اثر مهم‌ترین عوامل موثر بر رفتار سینتیکی واکنش نیتریلهکردن... 
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شکل ٩‏ نتایج آنالیز پارامتریک ۸) دماهای متفاوت درصد نیتروژن 


نتیجه گیری 
در مطالعه حاضر مدل مثر بر اساس روش گروهی بررسی 
داده ([01۲) توسعه داده شد تا سینتیک تشکیل نیترید 
سیلیسیم را پیش‌بینی کند. چندین پارامتر موثر مثل زمان. دما؛ 
درصد نیتروژن اندازهُ نمونه و اندازة ذرات سیلیسیم به‌عنوان 
متغیر اولیه برای مدل‌سازی درصد تبدیل لحاظ شدند. صحت 
مدل ۳۳1 توسعه‌یافته از طریق پارامترهای خطای آماری 
ارزیابی شد. نتایج آنالیز عملکرد نشان می‌دهند که مدل 
توسعه‌یافته با 8/۹7۲ و 2-/٩۳‏ دقت قابل توجهی در 
تخمین درصد نیتریده‌کردن داراست. به‌منظور اطمینان از دقت 
مدل توسعه‌یافته. نتایج مدل 0۷1۳11 با نتایج آزمایش واقعی 
نیتریده‌کردن انجام‌شده در ۱۵۷۳ کلوین و همچنین دو مدل 
عمومی و شناخته‌شده در زمینة نیتریده‌کردن سیلیسیم شامل 
مدل لی (۱۹۹۷) و چنگ (۲۰۰۰) مقایسه شدند. نتایج 
نشان‌دهنده این واقعیت بودند که دو مدل لی و چنگ در 


مراحل اولیه واکنش تخمین‌های خوبی را از پیشرفت واکنش 


سال سی و دوم شماره دی ۱۶۰۰ 


احسان شاه‌محمدی- فرهاد گلستانی‌فرد- علیرضا می رحبیبی 


نشان دادند ولی با پیشرفت واکنش بعد از دقيقه ۷۵ 
تخمین‌های این دو مدل خطای قابل ملاحظه‌ای داشتند. 
به‌هرحال مدل 1۳۳1) تطابق خوبی با نتایج آزمایشگاهی در 
تمام طول انجام واکنش از خود نشان داد. به‌منظور تعیین 
مهم‌ترین پارامتر موثر بر درصد پیشروی واکنش. آنالیز گاما 


3 


نیتروژن مهم‌ترین پارامترهای مژثر بر سینتیک واکنش بوده و 
دماء اندازهُ ذرات سیلیسیم و اندازه نمونه در جایگاه‌های بعدی 
قرار دارند. به‌منظور اطمینان از هم‌راستایی مدل با واقعیت‌های 
فیزیکی, آنالیز پارامتریک انجام شد. نتایج به‌دست‌آمده 
صحت مدل 0۷۳ را در تبعیت از الگوهای فیزیکی تأیید 


صورت گرفت. مطابق نتایج آنالیز حساسیت. زمان و غلظت کردند. 


هت 
منامزهاممتمط) نامع امملمعمع<-عسنهتع مه ۰۳۲ و.ظ .۲ ومب۷۱۱8 قمع رم روهظ و .۲ بنااننن .۲ ب,اممانعز۷ 


2018(۰) ,002979-2987 و(5) ,47 ۷۵۱ معاه 0126۵۲۵0۱۱6۵۵ ۳۱۱۵ ول , 1۱006 ]۱۷۲ )1و ۸۵/59/4۳۲۲ 0۴ 
6 صالصه] 1 ]۰ رب با رفع008۵621۷ 200 وبا ب۱(00 ۷ ,۱۷۲ ب911۷72 ویک مه8تا90 و.۲ .1 ,06170 و.۲ ۷۲۰ .۲ ,۷16112 
هل ,اهنا اتعتادظ رم51/513/۲10۵,:/۳۵/۲۵۵0 60ل2مصصه متاجوم‌جصم 1 متا 10 ععه1 طمروعط0 خر 

۳۶۱26۲۵۲۱6۱۷۵۵۵5, ۷۵۱, 45, ۵۰ 2, 00. 910-916, )2016(۰ 

۵ 111008 91-۳۴00 طا ممتفعتصصط عون ۵۶ ممتاهع۱۷۵0111. ریق رتتامقاه 220 وم ب0۵0160۷۵مطک و1 بقاف 010[ 
7 ۱۲6۰ ,47 ۷۵۱۰ ,داه۱/۵2 0۴1۳۵61۲۵۳۸۱6 ۳۵9۱۱۵ ,۱86۲۵۷ اطونا طمت)ه)۲۵1 مه «اعطصمند‌اماو فاوعه ۷ مار 
(2018) ,3927-3933 00۰ 

ممنصصهمن م0 مممنازگ مامزم۲ ۵ ولومطط ره ی ر0ظ 200 و.ت) ,120 پیک بتالناً ,۷۷۰ بتانا رل ۷208۵4 و.۷ با ] 
2010(۰) ,998-1000 .00 ,9 ۱۲0۵۰ ,25 ۷۵۱۰ بوهعکعع۲0 ۱۷۵۱۱۵۵۵8 ۵۱۱0 ۸/۵۸۵6 ۲1۵09۲۲۵ 

عصصنطمه]۱ ۵۶ مناد مهن( رت رصع ۷ 2020 ,.۷۷ بتالنا وباا2 رع1(00 ویک ۷۷ بش .6۷ مطفلن لاک ,2۳08 
۱ ر؟6956عع 0 ۹۵۱۱۵/۵0۱ ۵0 ۸۱۵/۵۳۲۵6 ۱۱10106 ممع‌نانه مصنل‌صتین ۷۷۵۲ 1۱8۸0۱660 عون مد وعتافتتمامهتق) 
2014(۰) ,1477-1482 .00 ,11-12 0۰ ,29 

ص کممصصمما۷۵عن1 مه طمتهعوم متمصهرمن [ - ]۱۱۲ یناه ۱۷۲ به‌تتاتداون ۷ 20 ری ,محصمتتا رک 501 
1991(۰) ,741-744 .00 4 ۱۱0۵۰ ,6 ۷۵۱ ,عجعع۵ظ و۸۵ له وا۵ ۷/۵/۵۲ ,1۳ ,۷۵۱ رطه2[ 

و6 ۱۱0۰ م34 ۷۵ را ۵۵6۵ 0 09۵ 16 ,۵۵2۵10085 عاهاه -90110 10۲ ۱۷۲۵061 متامصنک و یک رتماتهن) 
(1961) ,2010-2015 00۰ 

٩11160 6۲‏ 0۶ ووععم۲ عصنلجهها ممتامععع 10۲ عمج متامصل! ظ. .۲ یاه1۳008 200 ر.م) .ظ باصا ,.ءظ۷۷ ول 
997(۰) ,1119-1131 .00 ,9 ۱۲0۰ ,17 ۷۵۱۰ ,0628 6۲۵۳۱۵۵ ۱۱۵۵۵0 ۵۴۱۱۱۵ م0 0۵20 

امه ۵96۲ ممعنلند ۵ معتاممت! ممتلقلتتط م1 وبا .۲ تقو 20 وبلط ,1-۳1 بتاملنا .۷۷ ۲۰ رعصفطل) 


171267711061:101160 ۸616, ۷۵۱, 354, ۱۲0۰ 1-2, 00. 71-80, )2000(۰ 


٩1, 9۱/۵6۵ ۹661۱66, ۷۵01 494, ۱۵۰ 2, 00. 83-94, )2001(۰‏ 0۶ طمتاه 1 له ام ر.ل روطه ۷ ر.ت) بلق2 ر.کر رع2 ۷۷ . 


6۵ ۸۳/۱۵۹ ,ده انامه ممومتانه رها رممتانه ۵۶ ماهلا ویو رهته‌طاناعصتطه 200 و.] پتا2تصطلااه بل بقتهآ۲ . 


9661166, ۷۵۱۰ 363, 00. 355-559, )2016( 


۵ 000 00۷۷۵۲ 11608رو 0۶ و۷10۲ هه مم26 نت وبا مالک 200 ویک 11۲02 و.ما و0ظ2ولاکل و۷ بتامتا2 وال م۳۷ .12 


.(2011) ,291-293 00۰ ,119 ,۷۵۱ ,ال 0 ۹۵0۵6 ۵۲۵۱۱۱۵ 0۴۱۱۵ هل ,00069 طاتقه ۲2۲ دنامتنه۷ 


سال سی و دوم شماره دو ۳۹2 


3 مدل‌سازی اثر مهم‌ترین عوامل موثر بر رفتار سیتتیکی واکنش نیتریده‌کردن... 


مفصما ۱۷۲۱0۵9 زفصصماوزی 20 وعع۱۳۵۲1 رواه۱۷]2)۵۲1 متمز]۱۷ م2 مصول :11 مصصتتان ۲۷ .له ام ویو مطالا وبا ,۷۷۵208 و.۷ ون 
۰ ۷0۲1 ۱۱6۷ ,رد۸ , 046زا جمعنار فنامطام:محصه طز ۸04060 مممئلته گه ممتامتطداط جم ماگ :صمتاهتوماه1 
(2011) ,435-438 


-صمنامعع؟ مامت ۵ وا تالجم اقصمط]" 0صح عصلوومع۲۳۲۵ ر.2 رک6ک۵] 2۳00 ویک وتا و.۲۷ بتامط2 و26 بلا2 . 


,066 ۵۳۵6 ۱۵۲۵۵ ۱۱۵ 0۴ ۱۵ص و2۲0۵ ۳0۱۳۷0۵۲ ز ۵۴ ۳۲۲۵۵۲ :1 .۱0106 طمعت1ز ۳08060 


۲۷۵, 89, 0۰ 11, ۵0. 3331-3339, )2006(. 


٩ 00۷/0۵۲, ۵‏ 0 ط0تا2تتالظ مطا 10 ووممم م۱011 ۸ یگ رعتقک 220 ,۳۱-۷۷ رمما و.ل کته مب ملک . 


2005(۰) ,85-91 .00 ,1-2 0۵۰ ,406 ۷۵۱ بگ ۱8۵۵۵۵ 00 66166 


-0110 10 ۱0061 عمتامتطه مرمع-0مامههتص لمح عم ممتامهع 20۵۵ 0 ممفه۵ ۵ و2۱۷ بط ۷۷ 200 ,.۷۲ بع0طو1 . 


(۵0.1031-1041,))1 ,7 ۱6۰ ,26 ۷۵۱ ,6206 8۵۱۵68 ۱۵۱۱۵۵ ۷۹19 امصه امحطعط)180 آومم۲6۵011 فهع 


4۱ ۲0۰ و31 ۷۵۱ ,0۵97۵ ۸۵۱۳ و قطمتامهع 0110و-عهع طا معاوماه و۵۲ 0مباطاتتونل 0۶ فلورزلفمم رک .5 بقتاقطاظ . 


,(1985) ,642-648 .م2 


تاه 0۲۵ع-0م)مهم۵ا قح ام0م ممتامههعمم2۵ 0۴ ط0فلهمطهن وب[ رتهتله 200 و۷ ۷۰ ) م۷۷ .۷۲ ب02طو1 . 


,7 ۲6۰ ,26 ۷۵۱ ,۱266 ۱۱8۵۱۱۵21۳08۵ ۱۵۱۵۵ ر قلو مهد تمحصعطاهمو-2۵8 ۲6۵20110828-1 ٩011029‏ ص1 2۵061 


00۰ 1031-1048, )1971( 


66۵ ,۰ 009۵027۷ عط۱۷1 0 ۵ ۱۷1 دظ۲۵2۵)10 50110--ووع 10۲ ۵001 لمتاامنااو ۸ .۷۷ 1۰ رفصه ۲۷ 290 رل رزاا۹26 . 


۳۲۱۱8۵1۱۱66118 9621066, ۷۵۱, 25, ۲0۰ 6, 00. 1091-1107, )1970(۰ 

۱۵۲۵0۵۵0 متیهمنمد۱۷۲۱۵0 ما ۵ ممتهانمته مانام‌مهمن وبا بط رفعصتصصوژ مه رگ .۵ رتمععن و۷۷۰ بتاک 
۰ 2 ۱۰ ,73 ۷۵۱ 6۵6 6۳۵۱۱6 ۸۱۱۱۵۱6۵0 1۱۵ ]0 ۵9۱۱۵ , ومتصصحتعن تا ممعنازد 60ع۲ع)ص1-صمتامه۳۴ 
1990۰ ,286-296 

۵ ,۹111600 0۴ 020100اته مط) مد م۳070 قصه مد ۵۶ معصمتاط1 مط 1 وبا .۲ راهان 2240 وبا ,16۲۷1۵06820۷16 
(1979) ,1265-1268 00۰ ,5 6۰ ,14 ۷۵1 ,9667166 01۵۵۲1۵15 

۲۵۵ ۲۳۱ 1 1۵2۲2020100 م0مطامن ۳۵ ۵۶ مصتام۷]۵0 ۱۱600۲1 ماع ۰۸1110121 رل ب۱۷۵161 200 ر.ظ بتکلم[1۷]2 
2016(۰) ,3822-3834 .00 و8 ۲0۰ ,45 ۷۵ وعاهز۱/۵۱۵۲ 0۵6۲۵۳۵6 ۵9۲۵ , فلع 

00800 مطا ۰۱8۵۵0۲۵0۲2 وبا ب۸۵ بتصه۷ ۲ 220 ریک ,۸۷۲ بتماماهتا0ا۸ ملظ .9 بطله90۵۷60 مبلط ,۱۷۵1260862۳۵20 
۲ ۲۱0۱ عتناه۵ م6۳ ] -طم۲ مطا عصتام۲۲۵0۲ 10۲ طمهمتممد امممتامانام‌صومن ]901 ه 0ج ۱۷۱۵۵61 0)۷۵تافجمن) 
6 0 0۱/۲۲۱۵ ,۱۵۷۵۱۵0۴۵۵۵۲ ۱۷۱۵۵ عطلووع۱۵6ظ ۶۵۲ تناو وتاجدمصصمن) هر :امالخ ععلامو فصو ۵۶ 
2017(۰) ,467-477 .00 ,1 .۲0 ,46 ۷۵۱ ,کا۵ز۵۵۲ 

5 0660۵۵۵0۵۵ 0۶۴ معا رویط وملانک 20 ویک وفتالقکا بیط و96۷6 .۲ م001061) .1۷ و2عاط .لا ,126160251 
,261۲۵۱۵۱۷۵۱۵۲۱۵۱ 0۴ 09۳۱۱۵ ,۱۷۱۵۱۵۵05۳ 1همعصصلا عصله ۲۲ ۱۷۲۵۱۵121 0016ع1عمصصعط 1 ۵ عمصممصز0] ۲۵ و۵ 
.(2015) ,2068-2073 .00 ,6 ۲0۰ ,44 

تامهم ۱۵0011 لهتنمه اهلمزناته ری ب) رعما 200 وبلط .۲ متطلک ون 2 ,۲۱۵ ب.ظ .ظ بتطهتع ت۲2 یو .۲ ,۴600 
0 رخ 204 0 ۵۶ اتازصاهاه عطا ما مساجتهم‌جصها اصمصصاجعا احعط مه عاحمصوهاه هطاوماله ۵۶ عمصهازممر هام۲ عط 


2015(۰) و83 00.175-1 ,107 ۷۵۱۰ ,56۵066 ۷۵۵۲۱۵5 0۳۱۵۵۵۱۵0۲۵1 2۵110۷5 حصاتصها1) 10 


نشریة مهندسی متالورژی و مواد سال سی و دوم شماره دو 


13. 


20. 


21۰ 


22. 


23۰ 


24. 


25. 


۳3 


احسان شاه‌محمدی- فرهاد گلستانی‌فرد- علیرضا می رحبیبی 


۵ ,۷۵۱9 ,کتک ۵1 1۳۵۳۹۵۵۵6 1۳ , فصماونگ امن ۵۶ 0۲۷عظ۲ امتططمطرا۳۵ظ. ریت ظ رمکلمم‌مطله۷] 
1971(۰) ,364-378 .00 ,4 ۱0۰ ,)۹ ,۷۵۱ ,وراه 

۶ 0۵۷۵100 تمطاتن ۵ فماط9ا۲(۳0 ربا .0 بظ رماممصطل! ۲۷2 ره .0۰ 200 ربق .0 رمصمصط ۲۷2 ریت) بظ رملممصط۷۵] 
2 0۰ ,10 ۷۵۱ ,۸۵ ۷۸۵ ۵6۵8۳۱ ۴۸۲۲۳/۵۸۷ ,۲۰ )ت2ظ .فحصطاتعهع۱ه عصذاصمط 0202 ۵۶ 0مطاعحظ حنامع فطا 
2000(۰) ,00۰187-194 

۶ ۳۵0۱000 0۲ فممتامنوظ ها ۵ ممتاه ۵( رب .و رلته۳۱۵۲2۵-120 220 و.؟ ر900001ط2ظ ریم بطه۵۷ک1 
,240 ۷۵۱۰ ,۱۱۱668 0۳۱۱9۵۵ ۵ ال رصطانمعله ۷۲5 مصنعنا فتمامصهنع؟ ممتام]۷(-صنامیت لومتمصزرظ 
2016(۰) ,910-930 .00 ,6 0۰ 

986(۰) ,269-284 00۰ ,100008 مفصتلامی و عصزعمصعصظ [ابت صا دمتافتاهاو مه هتانماهط۳:۵۱0 .لا .۷ .0 مطاتص 
۵ 0 مج ممتامنلصملن )صمتاهط وی رلع00مطفطلقظ 220 وی .9 ۲۱2267120 بظ بطه۵۲ک 
77,۰ ۷۵1۰ ۱۱۷۱۱۲۵۱۵۱۱۵ ۱2 وه چوما0ع 0۴۱8۵۳۱۵۵8 ۱0/۱21۱ 0212 ۳1 وه 0عوده )مرمحطرووموعه لتاصعا20 
2018(۰) ,849-865 .00 ,2 

,۱۷/۱2/۵۲۵۲ ۵-۲20۵ رد0 5۷۲۵۲ 011۱16۲۵۵۵ .0 ۰ ,۲00۵120 200 ,.) .۸۵ الاک و.لا روتا>21 17 
2012(۰) ,00./97-801 ,7 ۲0۰ ,27 ۷۵۱ روععجعع0ظ م8 م۷۵ که داهز۸۵۱۵۲ 

ک۵۵۵ 0 ۱۵و , 02000تانه مممناته ۶ دلورلممه متامصل! ۷۵ تاماتاممنای ره رقمطته ۷ 200 ریت یک رطمعع1ظ 
1993(۰) ,00.2999-3013 ,11 ۱0۰ ,28 ۷۵1۰ ,661166 

۱۵20010۵40 0 متنام2ناصهه معط ما مامعره ۱۲2۳90۲۲ وباط .ر م۱۷۵۲ 220 بت ۷۵۵۲۵2 ریک .0 مفعطع۲۱ 
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سال سی و دوم شماره دی ۱۶۰۰ نشریة مهندسی متالورژی و مواد 


3 مدل‌سازی اثر مهم‌ترین عوامل موثر بر رفتار سینتیکی واکنش نیتریده‌کردن.. 


نشریة مهندسی متالورژی و مواد سال سی و دوم شماره دو, ۱۶۰۰ 


